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Abstract 

We take advantage of peculiar properties of three dimensional incompressible tur¬ 
bulence to introduce a nonstandard Exact Renormalization Group method. A Galilean 
invariance preserving regularizing procedure is utilized and a field truncation is adopted 
to test the method. Results are encouraging: the energy spectrum E{k) in the inertial 
range scales with exponent —1.666±0.001 and the Kolmogorov constant Ck, computed 
for several (realistic) shapes of the stirring force correlator, agrees with experimental 
data. 


1 Introduction 

Exact Renormalization Group (ERG) [2]p|0 represents a powerful tool for a nonperturba- 
tive analysis of quantum/statistical systems. It allows us to obtain, via an evolution equation, 
an effective action Fa in which only statistical fluctuations of scales smaller than 1/A are 
computed. Solving this evolution equation for A —> 0 leads us to the effective action F, the 
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Legendre transform of the generating functional of connected diagrams |]l|. This quite new 
tool has been applied to QFT and critical phenomena, mainly in conjunction with a 
Derivative Expansion |^, that is a truncation of Fa up to a maximum number of derivatives. 

In this letter we present a modihed ERG method for 3 — d turbulence and we explore the 
consequences of a simple held truncation in solving the set of coupled differential nonlinear 
equations of the ERG how. The truncation takes into account the ehects of the turbolent 
dihusivity corrections as well as noise corrections due to the energy cascade. 


2 The model 

Let us consider an isotropic homogeneous model of incompressible turbulence for a Newto¬ 
nian huid of viscosity u and density p, where the ehects of boundary/initial conditions are 
taken into account by the introduction of a stochastic noise term fa{,x,t). The “stirred” 
Navier-Stokes equations for the velocity Va{x,t) and pressure p{x,t) are 

VtVa{x, t) - uV‘^Va{x, t) + ^daP{x, t) = /„(x, t) , t) = 0 (2.1) 

where the convective derivative Vt = {dt + vpdjs) introduces the nonlinear term which gen¬ 
erates the turbolent behavior. These equations can be simplihed by the introduction of 
a transverse projector Papi'V) = {dap — and by the transverse convective derivative 

Dt = dt + P(V) {V0d0) : 

DtVa{x, t) - z/V^na(f, t) = fa{x, t) . (2.2) 

The stirring force / is chosen be a Gaussian stochastic white noise held which is completely 
determined by its two order correlator. Since the stirring force should mimic the instabilities 
occurring near the boundaries, its correlation length must be of the same order of the size 
of the system. If L is such a typical length, the correlation function 

< >= Po,0{^F{q){2'KY5^^\q + q)5{uj + uj') (2.3) 

must vanish for q » 1/L and for q « 1/L. The precise shape of the spectrum F{q) is 
not much important, nevertheless a normalization condition due to the energy conservation 
must be satished |^ < T >= / -^^F{q), where S is the rate of energy dissipated by a unit 
mass of huid. 
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To handle the statistical properties in a easy way and to throw a connection with the 
nsual QFT techniqnes, we need an action S for snch model. Snch action is given by the 
Martin, Siggia and Rose fnnctional [0], written in terms of the transverse held v and its 
canonical conjngated it: 

S'[h, n] = / d^xdtita{x,t)DtVa{x,t) — ita{x,t)h''V'^Va{x,t) 

+ d?xd^ydtUa{x,t)Paf}{V)F{x - y)ui 3 {y,t). (2.4) 

The generating fnnctional of connected Green fnnction can be written as 

W[J, J] = —i log ViiDv exp (iS[u, v] — iJv — iJu^^ . (2.5) 

where Ja{x,t) and Ja{x,t) are two anxiliary sonrces. We will compnte W as the Legendre 
transform of the effective action r['u,n], that is 

W[J, J] = r[h, v]-Jv- Jit (2.6) 

where it and v are determined by Ua{x,t) = and Va{x,t) = The effective 

action r[h, n] will be obtained by means of the ERG evolntion eqnation. 


3 ERG in turbulence 

In this section we will recall the main featnres on ERG method and we will apply it to the 
simple model of the previous section. The ERG prescription is based on the following points. 

• Gonsider a quantum/statistical theory described in terms of some held 0 and action 
S'[0]. For notational simplicity let 0 be a scalar held. If we modify the system in such 
a way that no quantum/statistical huctuations are allowed at distances greater than a 
certain scale, say, I/Aq with Aq very high, than the total contribution of the allowed 
huctuations is very small (and goes to zero into the limit Aq —> cxd). This means that, 
for this modihed dynamics which is described by some new action S'aq = S + AS\q, 
the ehective action Faq computed with Raq slightly dihers from Raq itself. 

• Given AS\ one can hnd the evolution equation for the ehective action Fa with respect 
to the scale parameter A. Let 

Wa[J] = -i log (^J V(j) exp (iS'A[0] - iJ4>)j (3.7) 
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be the generating functional of the connected Green function of this modihed theory. 
Following Wetterich [|^ we have: 

A^aFaI^ = A^aWaI j = A^a < A5a[0] > (3.8) 

where ^{x,t) = is the mean held of the modihed system. We can introduce 

also rA)^], the ehective action which collapses to S if huctuations were removed: 
rA[*h] = fA)^]—AS'A[‘h]. If AS'a[ 0] is a quadratic term AAa)*/*] = | / df‘xd'^y(j){x)Rf^{x— 
y)(j){y), equation (|3.8| ) becomes |Q 

ASaFaI, = hr (aSa-Ra ■ f wA-ir)') ■ (3-9) 

This is the evolution equation we are faced to solve. 

• Let us now consider the initial and boundary conditions for Fa. If we start with a 
sufficiently high cutoh Aq, huctuations are strongly suppressed and we can safely say 

if A = Ao , then Fao [<h] = S'[<F]. 

Since our goal is the computation of the “physical” ehective action F[<F], that is the 
sum of all the 1-particle irreducible graphs, we require that 


RmA^orA[<&] ^ r[<i>]. 


(3.10) 


Note that when lim^^oAS^ 
such that 


0 equation ( p.l0|) is satished. 


To sum up, if AS'a is 


1) for A —>■ Aq no quantum/statistical huctuation occurs and 

2 ) limA^oASA = 0 , 


then the solution of the evolution equation (|3.9|) gives us the ehective action of the 
theory. 


Let us now apply the method to the isotropic model of turbulence. We hrstly must look 
for AS A terms to add to the MSR action in order to satisfy the boundary conditions listed 
above. To do this, we need to discuss some physical aspects underlying the turbolent system. 
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The dynamics is characterized by the presence of two typical scales. One of them cor¬ 
responds to the macroscopic size L of the system. L is the scale at which the energy is 
introdnced into the system. Snch energy is transferred to smaller scales via the energy cas¬ 
cade driven by the nonlinear terms [Q [Q . This energy transfer stops at very smaller scales. 
The second intrinsic scale is the typical scale 7]^ at which snch energy is dissipated by the 
viscons term. The smaller scale rjd is called the “internal” or the “Kolmogorov” scale. By 
dimensional argnments it is easy to show (see for example [§]) that the internal scale is 
related to the physical parameters u and S by 

/ 3 \ 1/4 

Vd ~ ~ 

where TZ is the Reynolds nnmber of the model IZ = f/L/z/, being U the typical velocity at 
the scale L. Since IZ is the typical ratio between the inertial term {UV)U and viscons one 
it measnres the “strength” of the nonlinearity occnrring at the scale L. 

Tnrbolent regime can be realized only if TZ is mnch greater than nnity. If the Reynolds 
nnmber is very high, the system develops a wide domain of the Fonrier space, called “inertial 
range”, in which nonlinear terms dominate. Clearly, for momenta (1/T) << q « (l/'z^d) ~ 
7^3/4( 1 / 2 ^) bonndary terms (the stirring force in onr model) and viscons terms can 

be neglected. In this snbrange of Fonrier space we expect the system to be nniversal, self 
similar, homogeneous in space and isotropic 0 . 

Statistical properties of the inertial range have been studied intensively with many meth¬ 
ods (for a review see, for example, |P]). A way to obtain important information on the 
system is to apply a dimensional analysis (Kolmogorov, 1941 [§]). A simple analysis shows, 
for example, that the energy spectrum E{k) = < Ua{k,t)ui 3 {—k,t) > must follow 

a power law of the form E{k) = CK£‘^^^k~^ (Kolmogorov law), where Ck is a constant 
(Kolmogorov constant) of the order of unity. Clearly, being such analysis essentially a mean 
field approximation, we expect the exponents of the various moments be corrected by some 
anomalous dimensions. Such corrections are linked to the intermittent behavior of tnrbolent 
systems, that is strong nonlinear, rare, events usually associated to instantons 0. 

Let us try to use such informations to construct a regularizing term AS\. First of all 
note that if a system with a Kolmogorov momentum scale Kd = 1/rjd is stirred at momentum 
scales higher than Kd, no energy cascade develops. This can simply be inferred by looking 
at the Reynolds number TZ ~ {L/r]dY^^ << 1. Such a system is then quite non fluctuating 
and the effective action doesn’t differ too much from the bare action. This gives us the 
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right initial condition mentioned in point 1). Therefore we can introduce a AS'a term which 
modify the characteristic size L of the system into rj = 1/A, smaller than the internal scale 
r]d. Since the size L characterizes the force-force correlation function only, the AS'a term 
can be chosen as a functional that replaces the stirring force centered at the scale L with a 
stirring force at the scale rj = (1/A) < rjd « L. Given F{q), a possible choice of AS'a is 
then 

AS'a[0] = ^ J (fxd^ydtUa{x,t)RA^ai3U/3{y,t) 

= ^ J d^xd^ydtUa{x,t)Pay{V) {Fj^{x - y) - F{x - y))ui 3 {yR ), (3.11) 

where F\{q) = 2DoA~^h{q/A) has been obtained from F{q) by a rescaling L ^ rj = 1/A. 
The parameter Dq is linked, via the energy conservation, to the rate of energy dissipation £. 
By taking into account equation the normalization equation we have, in fact 

J dqq^h{q/A)^ = ^£ . (3.12) 


The regularized MSR action of this system S'a = S' -|- AS'a is now 

S'A['U,n] = / d^xdtUa{x,t)DtVa{x,t) — Ua{,X,t)h''V‘^Va{x,t) 
+ ^Jd^xd^ydtUa{x,t)Pay{V)F/^{x - y)up{y,t) , 

and the evolution equation (|3.9|) reads 

ASaFa = yr (aSa^a ■ (j^) Ia.a) ■ 


(3.13) 


(3.14) 


In our notation \u,u nieans the full, helds dependent, regularized propagator of the u 

held. 


Let us now check the boundary condition mentioned in point 2). The regularizing term 
( p.llp satishes /zttia^oAS'a = 0 in the following sense. Firstly we take the limit A —> 1/L so 
that ASi/a = 0 identically. Secondly the inhnite volume limit L —>■ oo is considered. 


Since our choice of AS'a satishes both initial and boundary conditions, the ehective action 
computed by the ERG how will be, in the limit A —1/L —0, the “physical” ehective action 
of the model. 


It should be noted that our choice of the regularizing term AS'a is somewhat unusual. 
In all the previous applications of ERG the AS'a term is introduced in such a way that all 
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the propagators of the theory are damped for q < A. This is not the case in our approach, 
since the correlator of the regularized system < uv >a has support on the entire Fourier 
space. Note, however, that all the correlation functions of the physical fields, i.e. < vv >a 
are suppressed for q < A both at the bare and at the “dressed” level. Indeed, the correlator 
< vv >A~ is proportional to the noise correlation function rtiu(p) Ri (Fa(p) + Ga(p)), 
which contains the correction of the noise term G^. The key point here is that the noise 
G'a(p) correction vanishes for p << L, so that the “dressed” noise term is damped for q < A. 
This is a consequence of the structure of the nonlinear terms in 3 dimensions and is deeply 
linked to the features of the energy cascade [|^, where the energy is mainly transported 
from low to high momenta and not vice-versa ( 2-d turbulence and situations which give rise 
to backscattering p[ are here not taken into account). 


4 The ERG flow 


The evolution equation (|3.14|) brings to an inhnite system of coupled nonlinear differential 
equations, one equation for every vertex = ^5fU=o- Because of the symmetries and 
of the structure of the ERG flow, however, some restrictions are found. Clearly the MSR 
action of the systems is invariant under Galilean transformations of an inhnitesimal “boost” 


Ca- Since the regularizing action AS'a is Galilean invariant too |^, one hnds Ward-Takahashi 
Identities (WT), which must be satished by Ta: 


d^xdt < {tCydpVa{x, t) — Ca) 


6Va{x, t) 


+ tCpdyUaix, t) 


6Ua{x, t) 


rA[h,n] = 0. (4.15) 


For a more detailed analysis of WT identities see The structure of the evolution equation 
poses some important restrictions on the form of Fa[m, n]. Indeed, the presence of the AcIaRa 
term in the evolution equation ( p.l4|) implies that no vertex correction with v legs only can 
be generated by the ERG flow . 

Even if some restrictions on the form of Fa[m, n] are found, the system of ERG differential 
equations is still of inhnite order. In order to handle the ERG equations, in this work we will 
explore the consequences of a held truncation in Fa. The simplest truncation of Fa which 
preserves WT identities is the following 

FU-nnc. [^, .y] = j d^xd^ydt {6^^\x - y)ua{x,t)DtVa{y,t) - Ua{x,t)f{A,x - y)V‘^Va{y,t)} 

+ \ ^ d?xd?ydtUa{x,t)Pap{y) {Ff^{x-y) + g{A,x - y))up{y,t). (4.16) 
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This ansatz for Ta takes into account the noise correction g{A, x — y) (energy cascade) as well 
as the damping correction /(A, x — y). The evolution equation for reduces then to a 

system of two coupled nonlinear differential equations in terms of f{A,x — y) and g{A,x — y), 
with the initial conditions 


f(ka,x-y) = v g(ka,x-y) = 0 


(4.17) 


Let us start with the computation of the effective viscosity f{A,x — y), which is is dehned 


as 


/(A,p)/ = 




\u=v=0 


(4.18) 


6u{p, 0)6v{—p, 0) 

The computation of the truncated effective action (|4.16|) can be simplihed by a rescaling 
of the momenta appearing into the trace, that is the loop momentum q, and the external 
momentum p. We rescale /(A, q) too in the following way: 

^ = a ; y = l'^ fi^p) = 


By taking into account of the explicit and implicit dependence on the scale parameter A of 
both the f{A,p) and Fa{p) functions, the evolution equation for 0 a(?/) reads 

(a^a - ydy - I) 4>A{y) = ^ f dxf\ dt (1 - t^){tx^/y - 2txy - y‘^) (3 + xdy) h{x) x 

+ 2/^ + 2txy)(j)A{x) {4 >a {kjx"^ + 2/^ + 2txy^ (x^ + 2/^ + 2,txy) + 0A(a;)a^^)) (4.20) 

The computation of the noise correction is similar. The g{A, q) function is defined by 


£/(A,p) = 


^^rA[A,n] 

^ 6u{p, 0)6u{—p, 0) 


As before we rescale ^'(A, q) in the following way 


(4.21) 


g{A,p) = 2DoA ^xa^v) 


(4.22) 


so that the evolution equation for xa( 2/) reads 

{AOa - ydy - 3) xa ( 2 /) = ^ f dx f\ dt ((1 - F)y‘^{x^ + y^ + ‘2Fx^ + ?,txy) (3 + xdyy) h{x) 
{Ha [k/x^ + 2 /^ + 2fx2/) + Xa [kjx^ + 2 /^ + 2fx2/))) x 
(^(x^ + 2/^ + 2tx2/)^0A(a^)0A {k/x^ + 2/^ + 2.txy^ 

{^(pA {kjx^ + 2/^ + ‘2txy^ (x^ + 2/^ + 2txy) + 0A(a;)x^)) 
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As we have seen in section 3, the Kolmogorov scale Kd is related to the viscosity and to the 
rate of energy dissipation by the relation K^, ~ The viscosity v can be expressed, 

then, in terms of £ (and so Dq) and Kd a.s v k, . Remembering that the initial 

cutoff Aq must be chosen well beyond the dissipation scale Kd, we hnd the initial condition 
for the 0A and xk functions 



oo; XAo(l/)~0. 


(4.24) 


The system of differential equations has been solved numerically with initial conditions ( 4.24 ). 
Universality in the inertial range means that the exact shape of the noise correlator should 
not be much important for the inertial range statistics. We have verihed the validity of 
this claim, within our approximations, by solving the ERG flow with several different noise 
terms. Each noise is, of course, normalized following ( 3.12|) . In addition, we asked the noise 
to vanish for null momenta, in agreement with the A —1/L —0 prescription, and to 
decrease exponentially to zero for g >> 1/L. 


5 Results 


We report here the results of our numerical solution for the ERG flow. We started at 
A = Aq = 1 with initial conditions ( [4.24|) and solved the differential equations with an 
explicit hnite-difference scheme. The results are obtained after a large number of iterations 
(A << Aq). The numerical computations show that the solutions approach a hxed point 
0 a(i/) —^ 0(i/);Xa(i/) —^ xiy)- For ouch choice of the stirring term, this hxed point doesn’t 
depend on the exact value of the initial condition for the viscosity term. 


A possible form of the stirring term is h{pL) = exp {—KpY- In hgure (|l]) we report 
the results for the hxed point functions (j){y) and xiu) obtained with this noise. 

The functions (j){y) and xiv) ure computed in the range 10 ~'<y = A< 10^. Glearly the 
subrange y » 1 corresponds to the inertial range. As it was expected, in this region the 
adimensionalized ehective viscosity 0 and induced noise y functions reach a scale-invariant 
regime, so that we parametrize them as 


Y{y)y»i = ; x{y)y»i = ly , 


(5.25) 
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where a and j3 {a k, (5 k. +0.37) are anomalous exponents. The hrst check for the validity of 
our approach is the computation of the energy spectrum. The correlation function at equal 
long time t is computed as 


< Ua{p,t)u^{-p,t) > = 


Pa^{p) 


2Do(A ^h{p/A)+g{A,p)'j 
/(A,p)p2 


(5.26) 


By taking into account (|4.19|) and (|4.22D , in the inertial range region {y » 1) we obtain 


E{p) 


1 2 X(P/A) 

27r2^ 0(p/A)p2 





(5.27) 


so that Ck = ^ data of hgure ([|) we can obtain the plot of the 

Energy spectrum by a simple ratio between the y and 0 data normalized by ^ (^) ^ 
the plot d^) we have resumed the spectra obtained with several different noise shapes, while 
in the table T1 we list the values of Ck and of the slope S = —5/3 + {(5 — a) obtained with 
each noise. 


Noise 

S (±0.001) 

Ck (±0.002) 

Noise 

S (±0.001) 

Ck (±0.002) 

2 —x^ 

xe 

-1.666 

1.124 

^4 g-+ 

-1.666 

1.146 

9 — 

xe 

-1.667 

1.267 


-1.666 

1.660 

2 —T® 

xe 

-1.666 

1.417 

8 —X® 

xe 

-1.666 

1.767 

1 

1 

-1.667 

1.489 

^10 g-xio 

-1.666 

1.784 


-1.667 

1.624 

^12 

-1.666 

1.785 


Table Tl: Values of the slope S and of the Kolmogorov constant Ck for various stirring forces 


6 Comments 


To sum up, in order to apply the ERG method to turbulence, we proposed an unconventional, 
symmetry preserving, regularization procedure. In addition, we don’t need to use the noise 
correlator of the form F[q) ^ q~y which is used in quite all the previous RG approaches 
0000 an exception is given by Brax’s work |^). Such form is clearly unphysical. 


as it is widely stressed [|] JTS 


To test the method a held truncation is applied to the ERG how. The results obtained 
by this quite crude approximation are encouraging. As it is clear from the table and hgure 
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(§), the second order inertial range statistics shows the required scale-invariance with scale 
exponents compatible with —5/3 ( for the second order statistics the anomalous exponent, 
if any, must be very small ) as well as strong independence with respect to large scale 
energy input details. The values of the Ck constant approach the value Ck ~ 1-78 in the 
limit of very narrow noise form (for experimental results see for example 

This approach is worth to be further analyzed, for instance, with higher orders held 
truncations or, more likely, with nonpolinomial approximations schemes. 


In a next paper |]^ the complete analysis of such procedure will be reported and different 
approximation schemes will be considered. 
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